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We present simulation results for the 2-flavour Schwinger model with dynamical Ginsparg-Wilson 
fermions. Our Dirac operator is constructed by inserting an approximately chiral hypercube oper- 
ator into the overlap formula, which yields the overlap hypercube operator. Due to the similarity 
with the hypercubic kernel, a low polynomial of this kernel can be used as a numerically cheap 
way to evaluate the fermionic part of the Hybrid Monte Carlo force. We verify algorithmic re- 
quirements like area conservation and reversibility, and we discuss the viability of this approach in 
view of the acceptance rate. Next we confirm a high level of locality for this formulation. Finally 
we evaluate the chiral condensate at light fermion masses, based on the density of low lying Dirac 
eigenvalues in different topological sectors. The results represent one of the first measurements 
with dynamical overlap fermions, and they agree very well with analytic predictions. 
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1. Motivation 

In 1998 a neat possibility was found to preserve chiral symmetry (in a modified form) on 
the lattice ^ ^. It has been used extensively in quenched QCD, but due to its computational 
demands the applications of dynamical chiral fermions are still in an early stage.' Hence it is 
strongly motivated to develop suitable algorithmic tools, in order to arrive — within a few years — 
at results that can be confronted with the light hadron phenomenology. 

Present chiral QCD simulations are restricted to coarse lattices; in particular thermodynamic 
studies typically use Nt = 4 and a ^ 0.28 fm. On such lattices the standard overlap operator is 
non-local. Locality improves, however, if we replace the kernel by a truncated perfect hypercube 
operator In quenched QCD, the locality of the resulting overlap hypercube fermion (overlap- 
HF) operator persists on rather coarse lattices ||6|]. At present, dynamical hypercube fermion (HF) 
simulations for QCD are under investigation [Q]; they will indicate how far the above property 
still holds beyond the quenched approximation. Thermodynamic tests show already that the cutoff 
effects for the HF are pushed to high energy [^]. 

In this work, we explore the feasibility of dynamical overlap-HF simulations with an algorithm 
which is peculiar to this type of Ginsparg-Wilson fermions. Our testing ground is the Schwinger 
model (2d QED), a popular toy model with certain features similar to QCD.^ Qualitative differences 
from QCD are the super-renormalisability of the Schwinger model and the absence of spontaneous 
chiral symmetry breaking. We consider the case of two degenerate flavours. 

2. The overlap hypercube fermion 

The Ginsparg-Wilson Relation (GWR) is a criterion for a lattice modified, exact chiral sym- 
metry [^, which was discovered by studying the properties of perfect and classically perfect 
lattice fermions. Since those formulations involve couplings over an infinite range, a truncation is 
needed, which distorts the perfect symmetry and scaling properties to some extent. For the free, 



optimally local, perfect fermion [ IC] the truncation to a unit hypercube preserves excellent scaling 
[ jlTl ] and chirality It leads to the form Dhf,.vv = P^l{^ ~ y)yti + ^{^ ~ y)y ^ vector term plus 
a scalar term {x,y are lattice sites). \nd = 2 these terms involve only couplings to nearest neigh- 
bours and across the plaquette diagonals. We gauge Dhf by multiplying the compact link variables 
Ux,ii € ?/(l) along the shortest lattice paths connecting x and y (for the diagonal the two shortest 



paths are averaged) [12]. Thus we arrive at the operator Dhf.at(^^)> which describes the HF.^ 

Since Dhf is 75-Hermitian, Dj^p = 75Dhf75> the exact chirality (which got lost in the trunca- 
tion) can be restored by inserting Dhf into the overlap formula which yields the overlap-HF 
operator 

DovHfH= (l-^)Di + m, Di°ijp=l + 75-^, //hf = 75pHF - 1) • (2.1) 

Y^HF 

//hf is Hermitian and D^^^'^p fulfils the GWR in its simplest form, {D^°|jp, 75} = D^ovwy^^mnv ■ I" 
practice we evaluate this operator by means of Chebyshev polynomials — after projecting out the 



'Recent status reports on dynamical overlap fermion simulations in QCD are given for instance in Refs. 
^Earliest efforts to simulate the Schwinger model with dynamical overlap fermions were reported in Ref. 
^We are using here the HF version which is denoted as CO-HF (chirally optimised hypercube fermion) in the Ref. 



Ip^. This is optimal for our algorithm to be described in Section 3. 
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lowest two modes of //^p' which are treated separately. The polynomial approximation was driven 
to some absolute accuracy of e (see below), so we deal with DovHF.e- 

Compared to H. Neuberger's standard overlap operator Dn |@], we replace the Wilson kernel 
Dw by Dhf Since the latter is an approximate solution to the GWR already, its transition 
Dhf — > ^ovHF is only a modest chiral correction (in contrast to the transition Dw Dn). This 
property is illustrated in Fig. ^ which compares the spectra of Dhf and Dqvhf for a typical config- 
uration at m = 0.03 and j8 = 5 on a 16 x 16 lattice."* 

Due to its perfect action background, we 
expect for Dhf also a good approximate rota- 
tion symmetry and scaling behaviour, which 
is then likely to be inherited by Dqvhf thanks 
to the relation Dqvhf ~ ^hf- This relation 
further suggests a high level of locality for 
Dovhf> since it deviates only a little from the 
ultralocal operator Dhf- These properties have 
been confirmed before in a Schwinger model 
study with quenched configurations, where mea- 
surement entries were re-weighted with the 




1 1.5 2 

Re (Dirac eigenvalue) 



fermion determinant [12] 



Figure 1: The spectra for Dhf and for DqvHF (up to 
a gap) at m — 0.03, in a typical, dynamical configura- 
tion at)3=5onal6xl6 lattice. Since the spectra 
are similar, Dhf is a good approximation to Dqvhf. and 
therefore approximately chiral. 
In order to simulate such fermions dynamically, the standard Hybrid Monte Carlo (HMC) 
algorithm would use the fermionic force term 



3. A preconditioned 
Hybrid Monte Carlo force 



- „-l f/^-l ^QovHF , ^6ovHF ^-l 

V^^oyHF V^ovHF"^^^^ t tiovHF J ^ovHF r ) 



(3.1) 



where 2ovHF = Ts^ovHF is the Hermitian overlap-HF operator, and Ax_^ are the non-compact gauge 
Unk variables. However, this force term is computationally expensive, and in addition conceptually 
problematic due to the discontinuous sign function H^f / \J in 2ovHF> see eq. (2.1). 

We render the force term continuous and computationally cheap by inserting only approximate 
overlap operators in the term (|3.1|).^ For the external factors we apply an overlap-HF to a low 



precision e', and we use //hf instead of 2ovHF in the derivatives (although this could easily be 
extended to a low polynomial as well). 



ovHF,e' 



HF 



■ + ■ 



HF „-l 



ovHF,e' 3 4 ' ^ovHF.e' 



2ovHF,e' ¥ ' 



(3.2) 



The Metropolis accept/reject step is still performed with the high precision overlap operator DovHF,e- 
Hence the deviations in the force are corrected, and the only point to worry about is the acceptance 
rate. A simplification, which reduces 2ovHF.e' to /sDhf . was originally proposed in Refs. [13] 



^Throughout this work we use the Wilson plaquette gauge action. 

^Such a modified force might also be helpful to achieve topological transitions more frequently, but we have no data 



for comparison with the force (3.1). 
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which reported a decreasing acceptance rate for increasing volume (that work used the "SO-HF" 
of Ref. However, it turned out to be highly profitable — and still cheap — to conect the 

external factors to a low precision. We chose 

e' = 0.005 (force term) , e = 10"^^ (Metropohs step) , (3.3) 

which increases the acceptance rate by an order of magnitude compared to the use of Dhf tlirough- 
out the force term. Note that the force we obtain in this way is not based on a Hamiltonian dy- 
namics, but the way we deviate from it (by proceeding from /sDhf to 2ovHF,e') does manifestly 
maintain the area conservation. 

4. Results for the acceptance rate, reversibility, locality and chiral condensate 

We performed production runs on a 16 x 16 lattice at j3 = 5 with five masses: m = 0.03, 0.06, 



0.09, 0.12 and 0.24. We applied the Sexton- Weingarten integration scheme [14] with a partial 
(5t)^ error cancellation (where Si is the step size). The time scales for the fermionic vs. gauge 
force had the ratio 1:5, but we did not observe a high sensitivity to this ratio. Our statistics of well 
thermalised configurations, separated by 200 trajectories, is given in Table |l|. 



Since the force ( |3.2| ) tends to push the trajectory a bit off the hyper-surface of constant energy, 
we kept the trajectory length (between the Metropolis steps) short. We chose it as £ = 1 /8, which 
is divided into 20 steps (i.e. dz = 0.00625); this turned out to be a good compromise in view of the 
acceptance rate and the dynamics between the trajectory end-points. Fig. ^ shows the acceptance 
rate (on the left) as well as the total number of required conjugate gradient iterations per trajectory 
(on the right). As usual, heavier fermions are easier to simulate. However, even down to our hghtest 
mass of m = 0.003 we obtained a useful acceptance rate ^ 0.3. 



a, 

o 




0.25 




0.25 



Figure 2: The acceptance rate (on the left), and the number of Conjugate Gradient iterations per trajectory 
(on the right, including all operations) as a function of the fermion mass m, on a 16 x 16 lattice at j3 = 5 and 
trajectory length £ = 1/8 = 20 • 5t. 



To study the quality of reversibility, we moved forth and back with a variable number of steps, 
and measured the (absolute) shift of the gauge action, IA^gI- Fig. || shows our results for the 
precision of the reversibility, still at 5t = 0.00625, for the masses m = 0.03 and 0. 12 and 0.24. The 
level of reversibility seems satisfactory. As we increase the mass, it improves significantly only 
at m = 0.24, as we also observed for 5t = 0.005. Our cunent results do not hint at any positive 
Lyapunov exponent, though this cannot be considered conclusive yet. 
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We tested the locality in the usual way [p^], 
by applying Dqvhf on a unit source rjy and mea- 
suring the decay of the function 



/(r) = T {£»ovHF,.rv.r], 



2 

II 



Co 

<1 




(4.1) 

We first consider the free case and demonstrate 
that this decay is much faster for the overlap- 
HF operator than for the Neuberger operator, 
see Fig. Q (on the left). On the right we show 
that the decay is still exponential for our dy- 
namically generated configurations, which con- 
firms the locality (and therefore the sensibility) 
of our Dirac operator. In the range that we 
considered, the mass has practically no influ- 
ence on this decay rate. A previous quenched 
re-weighted study revealed that the overlap-HF 

operator has a much higher degree of locality than the standard overlap operator Dn [12]. 
is observed here as well, since Dovhf at j8 = 5 is still far more local than even the free Dn. We 
repeat that improved locality also holds for the overlap-HF in quenched QCD and it enables 
the installation of chiral fermions on coarser lattices than the use of Dn. 



number of steps of size Sr = 0.00625 

Figure 3: The reversibility precision with respect to 
the gauge action for a variable number of steps of 
length dt = 0.00625. We show the results for our 
lightest mass and the two heaviest masses. We do not 
see an obvious indication for a positive Lyapunov ex- 
ponent. 
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Figure 4: The locality of the overlap Dirac operators, tested by the decay of the function (4.1), against 
the taxi driver distance in lattice units. On the left we compare our overlap-HF operator to the standard 
Neuberger operator (with H — YsiDvj — 1)) in the free case. The plot on the right shows the exponential 
decay of {f{r)) based on our overlap-HF simulations with various fermion masses at )3 = 5. 

At last we address the chiral condensate; it has been studied in the 1-flavour Schwinger model 
with quenched configurations in Refs. [|l^. For our case of two degenerate flavours,^ analytic pre- 
dictions were obtained for m <^ 1/ at low energy [18, |T^. This is realised in our settings, 
perhaps up to the case m = 0.24. In particular, Ref. predicts L{m) = -(y/i//) ~ 0.388(/7i/j8)^/^ 
based on bosonisation, while Ref. [18] arrived at a slightly lower coefficient ^ 0.37. Both predic- 
tions are marked in Fig. § (below, on the right). 

^For related work in QCD with dynamical overlap fermions and 1 or 2 flavours, see Refs. [p^]. 
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m 


0.03 


0.06 


0.09 


0.12 


0.24 


# of confs. 


205 


235 


221 


458 


100 


# of confs. at V = 





79 





325 


1 


# of confs. at V = 1 


205 


156 


220 


133 


46 


# of top. transitions 





1 


2 


3 


5 


(■^v=o) 




0.129(3) 




0.108(2) 




(■^|v|=l) 


0.171(2) 


0.173(2) 


0.171(2) 


0.165(3) 


0.169(4) 






"•^^"-0.031 




"•^^^-0.011 





Table 1: An overview of our statistics at different masses and in different topological sectors. Below we 
display our results for the leading non-zero eigenvalue A of D^^J^^p (with jack-knife errors). The values of 
E(m) were obtained based on the ratio between (A) in the topological sectors with |v| = and 1. 
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m = 0.24 
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Figure 5: Above: histories (left) and cumulative densities (right) of the leading non-zero Dirac eigenvalue 
A at different masses and in different topological sectors. Below on the left: the chiral condensate E as a 
function of the ratio (A|v|=i )/(Av=o), according to Random Matrix Theory pi]]. Below on the right we show 
the E values that we measured at two masses, cf. Table and the predictions of Refs. ^ |l|]. 

Fig. H (above) shows examples for HMC histories (left, illustrating the level of de-correlation), 
and cumulative densities (right), of the leading non-zero Dirac eigenvalues A.^ For our evaluation 



of the chiral condensate we made use of a formula given in Ref. [ |21| ] (for the £-regime), which 
expresses L{m) as a function of the ratio between (A ) in the sectors with topological charge | V | = 
and 1. For the masses considered, these functions are plotted in Fig. ^ (below, left). We read off £ 
for the eigenvalue ratios that we measured at m = 0.06 and m = 0. 12, see Table [ij Fig. |5| (below, 
right) illustrates our results, which agree with the predictions of Ref. |T^, 19] within the errors. 



^We used the eigenvalues of D^^^^, stereographically projected on the imaginary axis, where we take the absolute 
value. This treatment worked well also in quenched QCD pXl^. 
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Conclusions: We tested a force preconditioned HMC algorithm for the simulation of dynam- 
ical overlap fermions. It is applicable for the overlap-HF, but not for the standard overlap fermion, 
since it is designed for the case that the overlap kernel is similar to the overlap operator. In the 
2-flavour Schwinger model we obtained a useful acceptance rate and a decent precision of the re- 
versibility. A high level of locality is confirmed. We measured the chiral condensate at two fermion 
masses, and we obtained values for L{m) consistent with analytic predictions at low energy. 

Since the way to evaluate £ presented here leads to relatively large errors — the slopes of the 
functions shown in Fig. ^ (below, on the left) tend to be steep — we are now going to consider 
different methods for this purpose, along with an enlarged statistics. 

We thank M. Hasenbusch for helpful advice, and S. Dilrr, A. Kennedy and J. Verbaarschot for comments. 
J.V. was supported by the "Deutsche Forschungsgemeinschaft" (DFG). The computations were performed 
on the IBM p690 clusters of the "Norddeutscher Verbund filr Hoch- und Hdchstleistungsrechnen " ( HLRN). 
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